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0^ ■ ABSTRACT 
Os ■ 

We consider the Comptonization of an isotropic radiation field by a thermal 
<^ . distribution of electrons with non-vanishing bulk velocity. We include all rela- 

. tivistic effects, including induced scattering and electron recoil, in the derivation 

^ . of a kinetic equation which is correct to 0{9^, (39^, (3'^9e), where (3 is the bulk 

^ . velocity (in units of c) and 9e is the ratio of the electron temperature to mass. 

. The result given here manifestly conserves photon number, and easily yields the 

CN ■ energy transfer rate between the radiation and electrons. We also confirm recent 

. calculations of the relativistic corrections to the thermal and kinematic Sunyaev- 

O , Zel'dovich effect. 

OO . 

Subject headings: radiation mechanisms: thermal — radiative transfer — scat- 

P ' tering 
^. 

• 1— I , 1. Introduction 

' Many astrophysical processes depend on the Compton scattering of photons by a thermal 

distribution of electrons. Examples from cosmology include the evolution of anisotropies in 
the cosmic microwave background (CMB) radiation in the pre-recombination era [(Peebles 



fc Yu 1970)1 , and the Sunyaev- Zel'dovich distortion [(Zel'dovich fc Sunyaev 1969)| in the 



CMB spectrum due to passage through hot clusters. The Comptonization is described by 
the Boltzmann equation, which is an integro-differential equation for the evolution of the 
photon distribution function. For a detailed discussion of the Compton scattering kernel, 
see Kershaw, Prasad, & Beason (1986). In situations of astrophysical interest, it is often 
the case that the energy transfer during scattering is sufficiently small to allow one to make 
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a Fokker-Planck expansion of the scattering kernel. In this case, the Boltzmann equation 
may be replaced by a hierarchy of differential equations for the moments of the distribution 
function. 

The most complete analytic treatment to date of Comptonization in a moving media 
was given by Psaltis & Lamb (1997), who derived the first two moments of the photon kinetic 
equation for arbitrary (anisotropic) photon and electron distributions. This derivation was 
correct to first-order in h,ujjm(.(? and = kBTf./nif.c^^ where uo is the photon frequency, T^. 
is the electron temperature and mg the electron mass, and to second-order in the bulk 
velocity of the electrons in units of c. Psaltis & Lamb (1997) showed the importance of 
using the full relativistic cross section for Compton scattering in order to obtain the kinetic 
equation to the required accuracy. Although the generality of the results in Psaltis & Lamb 
(1997) ensures that they are sufficient to describe many astrophysical processes, there are 
situations of interest where the inclusion of higher-order relativistic corrections is necessary. 
One such example is the Sunyaev-Zel'dovich effect for hot clusters, where /csTe may be as 
high as 15 keV. In this example, the radiation is initially isotropic in the cosmic frame, and 
the optical depth to Compton scattering through the cluster is sufficiently small that the 
effect on the scattering of anisotropics induced by any peculiar velocity of the cluster may 
be neglected (the probability of a photon undergoing multiple scattering is very low). The 
effect of relativistic corrections on the thermal Sunyaev-Zel'dovich effect {(3 — 0) has been 
studied numerically (Rephaeh 1995; Rephaeh & Yankovitch 1997) and analytically (Stebbins 
1998; ChaUinor & Lasenby 1998; Itoh, Kohyama, & Nozawa 1998). Recently, these analyses 
have been extended to include the effects of any pecuhar velocity of the cluster (Nozawa, 
Itoh, & Kohyama 1998; Sazonov & Sunyacv 1998). It was shown that for typical values of 
the peculiar velocity j3 ^ 1/300, the relativistic corrections to the kinematic effect are ^ 8%, 
and arise from a term of 0{I39^ which is not included in the analysis of Psaltis & Lamb 
(1997). 

In the course of their derivation of the corrections to the kinematic Sunyaev-Zel'dovich 
effect, Nozawa et al. (1998) described a covariant method for performing a Fokker-Planck 
expansion of the photon kinetic equation for 7^ by extending the method used in Challinor 
& Lasenby (1998), where the electron bulk velocity was neglected {(3 = 0). This method 
includes all relativistic effects, including induced scattering and electron recoil in a unified 
manner. However, Nozawa et al. did not go on to evaluate the detailed form of the kinetic 
equation including all these effects; instead they gave only the resulting equation for a Planck 
distribution of isotropic radiation at temperature in the limit that <^Tg. This equation 
is sufficient to describe the dominant corrections to the Sunyaev-Zel'dovich effect (which was 
the purpose of their paper), but since the effects of recoil and induced scattering do not enter 
in this limit, their results are not sufficiently general to describe other interesting properties 
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of the Comptonization process, such as the energy transfer rate between the electrons and 
(hot) radiation. In fact, the dominant corrections to the Sunyaev-Zel'dovich effect have been 
derived independently by Sazonov & Sunyaev (1998) with recoil and induced scattering 
neglected from the start. Their calculation is much simpler than that in Nozawa et al. 
(1998), since they need only use the Thomson cross section. In the present paper, which 
is intended to complement the paper by Nozawa et al. (1998), we derive the detailed form 
of the photon kinetic equation describing Comptonization of an initially isotropic radiation 
field in moving media, in the limit of low optical depth. In this limit, the Compton scattering 
term in the Boltzmann equation may be evaluated for an isotropic distribution, and, since 
multiple scattering is very improbable, the effects of the velocity-induced anisotropies on 
subsequent scatterings can be safely ignored. We feel that this result, which is omitted 
from the Nozawa et al. (1998) analysis, could be valuable to the astrophysics community at 
large since Comptonization is central to many problems. Relativistic effects may be fully 
included by a systematic expansion in the parameters 9(, and 13. The new result given here 
is written in a form that manifestly conserves the number of photons, and allows a simple 
calculation of the energy transfer rate between the electrons and the radiation. In the limit 
of /3 = we recover the expression given in Challinor & Lasenby (1998). For /3 ^ 0, our 
result yields corrections at higher order than given elsewhere the hterature. We give the 
kinetic equation correct to 0{9^, ^6^, f3'^6e), before calculating the rate of energy transfer 
for a Planck distribution correct to 0{6^,6^, (36^, (36r), where 6r = kBTj./meC^ . We end with 
a brief discussion of the terms in the kinetic equation that are required to describe the 
Sunyaev-Zel'dovich effect, obtaining results in full agreement with those in Nozawa et al. 
(1998). We use units with = c = fc^ = 1 in the following, unless stated otherwise. 



2. Extending the Kompaneets equation 

We consider the Comptonization of an unpolarised, initially isotropic radiation field by 
a thermal distribution of electrons at temperature Tg which has bulk velocity f3 relative 
to the radiation. The optical depth is assumed to be sufficiently low that the radiation 
may be treated as isotropic throughout the interaction with the medium. We shall work 
exclusively in the frame in which the radiation is isotropic, but it should be emphasised 
that we express our final results in terms of the electron number density in the rest frame 
of the scattering medium, N^. The electron distribution function is denoted by f{E,p) and 
the photon distribution by n{u!, k). (Note that we use distribution functions normalised to 
equal the mode occupation numbers.) Here, the electron energy is E, and the direction of 
propagation p. For the photons, u denotes the energy and k the direction of propagation. 
Neglecting the effects of electron degeneracy, the Boltzmann equation for the evolution of 
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k) may be written as (Buchler fc Yeuh 1976) 



Dn{uj, k) 
Dt 



-2 



(2^ 



d^p'd^k'w[n{iu, k) [1 + n{iu', k)]f{E, p) 

-n{u;', k)[l + n{u;,k)]f{E',p')}, 



(2-1) 



where the operator D/Dt denotes dt + k-'V. The invariant transition amphtude W for 
Compton scattering of a photon of 4-momentum k^^ by an electron (of charge e and mass 
rUe) with 4-momentum p^, to a photon momentum k'^ and an electron momentum p'^^ (whose 
energy is E') is KBerestetskii, Lifshitz, fc Pitaevskii 1982) 
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with n = —2p^k^ and k' = 2p'^k'^. 

The electron distribution function is assumed to be a relativistic Fermi distribution in 
the frame moving at /3. Since we are ignoring degeneracy effects, in the frame of the radiation 
we have 

f{E,p) ^e-^^^^~f^-P^-f''^^/^% (2-4) 

where fXe is the electron chemical potential and 7 = (1 — Substituting for f{E,p) 

into equation (P^), setting n{u, k) = n{uj) in the integrand, and expanding the distribution 
functions in powers of Ax, where 

X = uj/T„ (2-5) 
Ax = {uj'-uj)/Te, (2-6) 

gives the Fokker-Planck expansion for an isotropic radiation field |(JNozawa et al. 1998 )| 
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We calculate the I^^i and coefficients by expanding the integrands of equations ( P-^ ) 
and ( |2-9| ) in powers of p/rrie and uj/rrie. These integrations are ideally suited to sym- 
bolic computer algebra packages (we used Maple). To derive a kinetic equation correct to 
0{9^, (39^, P^Oe), one must evaluate Ji,o through to Js^s, and Jo through to J4. Substituting 



the resulting coefficients back into equation 
1 Dn{uj,k) 
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7 d i 



dx I 
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we find the kinetic equation 
5 
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where /i is the cosine of the angle between the photon momentum and the peculiar velocity of 
the electron distribution, = k-^/P, the Pi{fi) are the Legendre polynomials, and is the 
number density of electrons in the frame where the bulk velocity vanishes. The coefficient 
Ci is given by 
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and C3 by 
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Equations ( |2-10| )-( ^-13| ) are the main result of this paper. Some of the terms in equa- 
tion (|2-10| ) have been given previously in the literature; setting (3 = recovers the kinetic 
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equation given in Challinor & Lasenby (1998) which we used to investigate corrections to the 
thermal Sunyaev-Zel'dovich effect (higher-order corrections to the (3 = equation were given 
by Itoh et ah (1998)), the 0{(3^) term inside the curly braces is implicit in the / = moment 
equation given in Psaltis & Lamb (1997) (but not the 0{j3'^6e) term which is fourth-order 
in the electron velocity), and the 0(/?^) term and part of the 0{(39e) term are implicit in 
the analysis of Sazonov & Sunyaev (1998), although they have ignored the parts of these 
terms arising from induced scattering and recoil effects. For any particular application, the 
validity of neglecting the higher-order terms in equation (|2-10|) should be carefully checked. 
The terms that we have given are more than sufficient to describe the kinematic Sunyaev- 
Zel'dovich effect for typical cluster parameters, (3 ~ 1/300 and ksT^ — lOkeV. (In Nozawa 
et al. (1998) it is shown that the 0(/3^) terms are insignificant for these parameters, while the 
0{(36e) term gives a correction of —8.2% to the kinematic Sunyaev-Zel'dovich effect at the 
position of the zero of the thermal effect, and the 0{f36^) term gives a correction of +1.3%.) 
If required, higher-order terms in the kinetic equation can be derived from the Fokker-Planck 
expansion (eq. (^^)), although in practice the evaluation of the I^^i and Jk rapidly becomes 
prohibitive. 

We have written equation ( |2-10|) in a form that manifestly preserves the total number 
of photons, as required for Compton scattering. For /? = 0, we obtain a generalisation 
of the diffusion approximation to the Boltzmann equation (see, for example, Prasad et al. 
(1988)). Note that we have derived our equation (eq. ( |2-10D ) by a systematic expansion of 
the original Boltzmann equation (eq. (|2-1|) ) in l/me, we have not appealed to the heuristic 
arguments that form the basis of the diffusion approximation. For a spatially localised, 
isotropic distribution of photons, the time rate of change of the total photon number A^^ is 
given by 



dt J (27r)3 Dt 

where d^x is the spatial measure, and the factor of two accounts for the two polarisations. 
Integrating equation ( |2-10|) over photon momenta, the terms involving /i vanish by virtue of 
the integral over solid angles, and the //-independent terms (those in curly braces) vanish after 
integration over photon energies, since these terms are written in the form of a conservation 
law. It is not hard to show that equation ( |2-10D with (3 = admits static, homogeneous 
solutions with n = l/(exp(a; — z/) — 1) as required for thermodynamic equilibrium. 
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3. Rate of energy transfer 



The rate of increase of energy density E,. in the radiation due to Compton scattering is 
given by 
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Substituting for Dn/Dt from equation ( |2-10|) , we see that only the /i-independent terms 
contribute to the energy transfer. Substituting a Planck distribution at temperature for 
n and performing the integral, we find 
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where ({x) is the Riemann Zeta function, and 6r = Tr/rrie. The terms in the first square 
bracket in equation ( p-2| ) are independent of /?; they tend to equalise the radiation and 
electron temperatures. These terms (which were also given in Challinor & Lasenby (1998)) 
were first derived by Woodward (1970), where higher-order terms were also given. The 



terms in the second square bracket in equation (|3-2D represent the lowest-order effects of 
the electron bulk velocity on the energy transfer. The first such term AErO'TNeiS'^ /3 is well 
known (see, for example, Sazonov & Sunyaev (1998) and references within). 



4. The Sunyaev-Zel'dovich effect 



For CMB photons passing through a cluster at redshift z, the average value of a; = u/Tg 
is X ~ 6.2 X 10~^{1 + z)/kBTe, where ksT^ is expressed in eV. Since the electron temperature 
for a hot cluster is typically ~ lOkeV, it follows that a; -C 1. In this limit, equation ( P-IQ ) 
reduces to 
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Setting n — l/(exp(Q;j:) — 1), where a = T^/Tj. recovers the combined thermal and kinematic 
Sunyaev-Zel'dovich spectral distortion given in Nozawa et al. (1998) and Sazonov & Sunyaev 
(1998), where the relative importance of the various terms for typical cluster parameters is 
discussed in detail. 



5. Conclusion 

Using the covariant Fokker-Planck expansion described in Nozawa et al. (1998), we have 
derived a kinetic equation describing the interaction of an isotropic radiation field with a 
thermal distribution of electrons, which moves at bulk velocity cfi relative to the radiation, in 
the limit of low optical depth. Relativistic effects are included to 0{6^, (36^, j3'^6e), or equiv- 
alently to 0(9^, /39^, /3^0) where 6 is either of hu/meC^ or kBTe/rrieC^, which is sufficient 
to describe the corrections to the kinematic Sunyaev-Zel'dovich effect for typical cluster pa- 
rameters (Nozawa et al. 1998) The method may be easily extended to include higher-order 
relativistic effects if required. We have calculated the rate of energy transfer between a 
Planckian radiation field and the electrons, obtaining the usual "thermal" and "kinematic" 
terms, as well 0{l3'^9) "interference" terms. Speciahsing to the hmit <C T^, we confirm 
the relativistic corrections to the thermal and kinematic Sunyaev-Zel'dovich effect given in 
Nozawa et al. (1998) and Sazonov & Sunyaev (1998). 

We would like to express our gratitude to Roberto TuroUa and Silvia Zane for bringing 
a number of useful references to our attention. 
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